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The STEREO/WAVES experiment (SWAVES) onboard the two STEREO spacecraft 
(Solar TErrestrial RElations Observatory) launched on Oct 25, 2006, is dedicated to the 
measurement of the radio spectrum at frequencies between a few kHz and 16 MHz. The 
SWAVES antenna system consists of 6 m long orthogonal monopoles designed to measure 
the electric component of the radio waves. With this configuration direction finding of 
radio sources and polarimetry (analysis of the polarization state) of incident radio waves is 
possible. For the evaluation of the SWAVES data the receiving properties of the antennas, 
distorted by the radiation coupling with the spacecraft body and other onboard devices, 
have to be known accurately. In the present context, these properties are described by the 
antenna effective length vectors. We present the results of an in-flight calibration of the 
SWAVES antennas using the observations of the non-thermal terrestrial Auroral 
Kilometric Radiation (AKR) during STEREO roll maneuvers in an early stage of the 
mission. A least squares method combined with a genetic algorithm was applied to find 
the effective length vectors of the STEREO-B/WAVES antennas in a quasi-static 
frequency range ( L antenna A waV e) which fit best to the model and observed AKR 
intensity profiles. The obtained results confirm the former SWAVES antenna analysis by 
rheometry and numerical simulations. A final set of antenna parameters are recommended 
as a basis for evaluations of the SWAVES data. 


1. Introduction 

The complex radiation coupling between the an¬ 
tennas aboard a spacecraft and other electric de¬ 
vices or metallic structures alters currents induced in 
the electric antennas by an incident electromagnetic 
wave, and, therefore, causes the distortion of the 
reception properties. As result the effective length 
vectors, which describe the main antenna reception 
properties such as directional characteristics and ef¬ 
fective length, differ from the expected ones based 
on properties of the ’’stand-alone” antennas. There¬ 
fore, for the accurate evaluation of the data acquired 
by the radio instruments onboard spacecraft the re¬ 
ception properties of the antennas, influenced by the 
spacecraft body, have to be known precisely enough. 
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This can be done by means of several well devel¬ 
oped techniques [ Macher , 2008], such as rheometry 
[Rucker et al., 1996; Macher et al., 2007], aneclioic 
chamber measurements [Riddle, 1976], in-flight cali¬ 
bration [Vogl et al., 2004; Panchenko , 2004; Cecconi 
and Zarka, 2005] or computer simulations [Oswald 
et al., 2009; Rucker et al, 2011]. 

The rheometry method uses the scaled model of 
the spacecraft suspended in an electrolytic tank and 
the antenna reception properties are determined by 
measuring a response of the electric antennas to the 
electric field in the quasi-static frequency range, i.e. 
when the radio wavelength is much greater than 
the antenna length (short electric antenna) [Rucker 
et al., 1996; Macher and Oswald, 2011]. A very pow¬ 
erful technique to analyze the reception properties of 
the spaceborne antennas is a computer-based wire or 
patch-grid simulation in which the spacecraft body 
and antennas are modeled as mesh of wires or patches 
[Fischer et al., 2001; Oswald et al., 2009; Sampl et al., 
2011, 2012]. Numerical solutions of the electric and 
magnetic field integral equations by means of numer¬ 
ical electromagnetic codes yield the satisfying deter¬ 
mination of the antenna effective length vectors and 
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antenna radiation pattern for quasi-static as well as 
for higher frequency range. These computer simula¬ 
tions allow to predict the reception properties of the 
antennas before launch of the spacecraft. 

The main idea of the in-flight calibration is the 
determination of the reception properties of the an¬ 
tennas after spacecraft launch using the natural radio 
sources such as galactic background or planetary ra¬ 
dio emission. An in-flight calibration procedure con¬ 
sist of two main parts: 1) determination of the an¬ 
tenna system gain including effective lengths of the 
antennas as well as base and antenna capacitances, 
and 2) determination of the antenna directivity - i.e. 
the directions of the antenna effective lengths vec¬ 
tors. 

The antenna gain and effective lengths can be 
determined using the quasi-isotropic non-thermal 
galactic background as a reference radio source. This 
method has been implemented to derive the abso¬ 
lute flux density measurements of the Cassini/RPWS 
[Zarka et al., 2004] and the effective length of the 
STEREO/WAVES electrical antenna system [Za¬ 
slavsky et al, 2011]. 

The effective antenna length vectors of the electric 
antennas can be investigated by analyzing the tem¬ 
poral variation of the intensity of the radio emission 
emitted from the point source with known location. 
In this method the spaceborne antenna system must 
rotate with respect to the direction of the incident 
radiation. The in-flight calibration method is one of 
the most reliable techniques to determine the recep¬ 
tion properties of the antenna because the method 
deals with the real observations. The methodology 
of this calibration is well described in Vogl et al. 
[2004]. The method has successfully been applied to 
derive the antenna properties of radio experiments 
onboard e.g. ISEE-3 [Fainberg et al., 1985], Voy¬ 
ager [Lecacheux and Ortega-Molina, 1987; Wang and 
Carr, 1994], Interball-2/Polrad [ Panchenko , 2004] or 
Cassini [Vogl et al., 2004; Cecconi and Zarka, 2005] 

In this paper we present the results of determin¬ 
ing the effective length vectors of the antennas of 
the WAVES instrument onboard STEREO space¬ 
craft. We used the terrestrial Auroral Kilometric Ra¬ 
diation (AKR) - intense radio emission from auroral 
regions [Gurnett, 1974], as a radio source. The an¬ 
tenna effective length vectors have been determined 
by fitting of the model-predicted temporal variations 
to the time profiles of the AKR intensity measured 
during STEREO/WAVES roll maneuvers. The ob¬ 


tained results are compared with antenna effective 
length vectors evaluated by methods of rheometry 
and numerical computer simulation. 

The STEREO/WAVES radio experiment is de¬ 
scribed in chapter 2. The method of in-flight cali¬ 
bration of the SWAVES antenna is presented in sec¬ 
tion 3, and the results and error analysis are given in 
section 4. In section 5 we discuss the results. 

2. STEREO/WAVES antenna system 

Solar TErrestrial RElations Observatory (STEREO) 
consists of two identical spacecraft (STEREO-A and 
STEREO-B), launched on October 25, 2006. The 
main goals of the STEREO mission are investigations 
of the three-dimensional structure and evolution of 
the solar coronal mass ejections (CMEs) as well as 
the study of the CMEs interaction with the Earth’s 
magnetosphere. The set of scientific experiments on¬ 
board STEREO provide complex measurements of 
the electromagnetic and local plasma waves, as well 
as 3D magnetic field components, and plasma pa¬ 
rameters like the solar wind bulk velocity, density 
and temperature. After several highly eccentric geo- 
centrical orbits and lunar swing-by maneuvers (Jan¬ 
uary 2007) the two STEREO spacecraft split up and 
were inserted into a heliocentric orbit. One of the 
spacecraft, STEREO-A (Stereo Ahead) leads Earth 
in an orbit which is slightly closer to the Sun, while 
the other spacecraft, STEREO-B (STEREO Behind) 
trails Earth orbiting slightly outside the terrestrial 
orbit. 

The SWAVES radio experiment is an interplane¬ 
tary radio burst tracker that observes the generation 
and propagation of the radio disturbances from the 
Sun to the orbit of Earth [Bougeret et al., 2008]. 
SWAVES consists of fixed, high (HFR) and low 
(LFR) frequency receivers as well as of a time domain 
sampler (TDS). The antenna system of the SWAVES 
instrument include three nearly mutually orthogo¬ 
nal electric monopoles X, Y and Z. Each of the 
monopoles has a length of 6 meters. 

The super-heterodyne swept frequency receiver 
HFR (High Frequency Receiver) includes two sep¬ 
arate parts (or ’’logical receivers”) HFR1 and HFR2 
which cover the frequency ranges 0.125 - 2 MHz and 
2-16 MHz, respectively, with 50 kHz frequency reso¬ 
lution. The time resolution depends on the work¬ 
ing mode. HFR1 can work in direction finding 
mode providing spectral and complex cross-spectral 
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power densities of incident waves. HFR.l is a dual 
channel receiver which can be connected to a pair 
of monopoles X, Y or Z or to combinations of 
monopoles and ’’pseudo dipoles” (monopoles X or 
Z can be paired with the Y antenna). Thanks to the 
two parallel analysis channels, HFR1 measures in¬ 
stantaneously four values: autocorrelation of each of 
the two antennas and the real and imaginary part of 
cross-correlation between pairs of the antennas (see 
Bougeret et al. [2008] for more details). 

In the so-called direction finding mode (DF) the 
HFR1 switches between the pairs of the antenna at 
each frequency step, providing quasi-instantaneous 
auto and complex cross-spectra density of the wave 
between different antenna combination. HFR1 is op¬ 
erated in two special direction findings modes: 

DF1 - ’’two antenna mode”, in which two possi¬ 
ble antenna configurations are used: 

DF1 ’13’ - one monopole Ex and one dipole 
Ezy (dipole consist of Ez and Ey monopoles). 
The output signals in this mode are: auto correla¬ 
tions {Ax, Azy) and cross-correlations (5f x/zy) and 
(^x/zy)- 

DF1 ’31 ’ - one monopole Ez and one dipole 
E\y (monopoles Ex and Ey are connected as a 
dipole). Operating in this mode HFR1 at each fre¬ 
quency provides four quasi-instantaneous output sig¬ 
nals: two auto correlations (Az, Axy ) and the real 
(^z/xy) and the imaginary (Qz/xy) part of the 
cross-correlation between Ez and Exy ■ 

DF2 - ’’three antenna mode”, uses sequence 
combination of all monopoles Ex, Ey, Ez- This 
mode yields nine quasi-instantaneous output signals: 
three auto correlations (Ay, Ay , Az)', three real 
parts (JRa'I’, ^yz, Kzx) and three imaginary parts 
(3xv, ^yz, ^zx) of cross-correlations. 

3. Method and data preparation 

3.1. AKR observation 

Several spacecraft roll maneuvers have been sched¬ 
uled at the beginning of the mission. During this 
maneuvers the spacecraft were at 90 — 140 Re 
away from the Earth. Each maneuver lasted 10 
hours and consisted of 10 consecutive rolls (6° per 
minute). In particular, STEREO-A performed three 
roll maneuvers: 1st on 12/18/2006 06:50-16:50, 2nd 
on 12/20/2006 07:20-17:20 and 3rd on 12/23/2006 
04:30-14:30. Second spacecraft, STEREO-B, per¬ 
formed only two roll maneuvers on 01/29/2007 05:10- 



Stereo B, 01/29/2007 (DOY= 29) 
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Figure 1. Dynamic spectra of the radio emission 
recorded by STEREO-B on 29 and 31 Jan., 2007 (pan¬ 
els a and c). Terrestrial AKR is observed in a frequency 
range up to 700-800 kHz. The panels b) and d) show the 
roll angle of the spacecraft in the GSE coordinate frame. 
This angle shows the times when STEREO-B performed 
roll maneuvers: on 29/01/2007 between 05:10-15:10 UT 
and 31/01/2007 between 04:50-14:50 UT. 


15:10 and 01/31/2007 04:50-14:50. These spacecraft 
roll maneuvers can be used to determinate the effec¬ 
tive length vector of each antenna. 

The in-flight calibration requires the observation 
of the radio emissions from a known position rel¬ 
ative to the rotating antenna system. During the 
roll maneuvers only one spacecraft, STEREO-B, ob¬ 
served the terrestrial Auroral Kilometric Radiation 
(AKR) with the signal/noise ratio sufficiently high 
to perform the antenna calibration. The reason why 
only STEREO-B observed AKR is that, after the last 
lunar swing-by maneuver, STEREO-B flew by the 
Earth crossing the dusk side of the magnetosphere 
(from dayside to nightside). Such a trajectory en¬ 
abled a long-lasting (several months) observation of 
the terrestrial AKR (see [Panchenko et al., 2009]), 
whose sources are known to be merely fixed in local 
magnetic time of the dusk-night part of the magne¬ 
tosphere. In the same time STEREO-A was on the 
dayside magnetic local time sector and, therefore, it 
was not able to observe the AKR. 

Fig. 1 shows the dynamic radio spectra recorded 
by STEREO-B during SWAVES antenna rolls on 
29 and 31 Jan., 2007. The bottom panels (b) 
and (d) show the roll angle of the spacecraft 
in the GSE (Geocentric Solar Ecliptic) coordinate 
frame. The roll angle as well as the spacecraft 
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attitude and position have been calculated using 
the NASA NAIF/SPICE toolkit - a system which 
provides spacecraft ancillary information (http : 
/ / naif.jpl.nasa.gov / naif /tool kit.html ). 

As shown in Figure 1 only in one episode, i.e. 
on 31 Jan. 2007, AKR was observed continuously 
during the whole roll maneuver, whereas on 29 Jan. 
2007 the radio emission was detected only during two 
short time spans, i.e. 06:30-08:00 and 14:15-15:20. 
Therefore we have decided to use only the observa¬ 
tions performed on 31 Jan., 2007 between 04:50 and 
14:50 UT. During these time intervals HFR1 was op¬ 
erated consecutively in DF1 T3’, DF1 ’31’ and DF2 
modes i.e. in each consecutive frequency sweep 
HFR1 was switched between DF1 ’13’, DF1 ’31’ and 
DF2 modes. 

The first step of the data preparation deals with 
a subtraction of the frequency dependent back¬ 
ground which consists of receiver and galactic back¬ 
ground noises. The level of the background has 
been determined in the same way as it was done 
for Cassini/RPWS calibration [ Vogl et al ., 2004, ap¬ 
pendix A ]. In particular, assuming that the back¬ 
ground is nearly constant with time and depends 
only on the frequency and working mode of the re¬ 
ceiver the background level was determined using a 
histogram of occurrence probabilities of intensities 
calculated for each frequency. The histograms were 
obtained during time intervals when no AKR activ¬ 
ity was detected. Then, the background level at each 
frequency was defined as the intensity which corre¬ 
sponds to the lower 5% occurrence level (left from the 
maximum of the occurrence histogram) (see Zarka 
et al. [2004]). The data points below the background 
were excluded from the analysis in order to avoid any 
bias in our data set. Since STEREO-B/WAVES was 
operated mainly in DF2 mode (in the beginning of 
the mission ) the frequency dependent background 
was calculated only for this mode. The background 
for the DF1 mode has been assumed to be at the 
same level as for DF2 mode, although, as was shown 
by Taubenschuss [2005] the background can be a 
slightly different for different antenna combination. 
It is worth to note, that the intensity of the AKR 
observed on 31 Jan. 2007 was 2 — 3 orders of mag¬ 
nitude higher than the averaged background level at 
frequencies below 1 MHz. Therefore, any potential 
inaccuracies in background subtraction do not play a 
significant role in the antenna calibration procedure. 


After the background subtraction, the intensity 
profiles have been established by an averaging of the 
data over the frequency bandwidth 225 - 475 kHz 
(6 frequency channels of the HFR1 receiver). These 
normalized intensity profiles are plotted as solid blue 
lines in Figure 3. The standard deviation of each 
mean value will be used as uncertainties of the ob¬ 
servations a.i in equation (3). 

3.2. Model predicted SWAVES outputs 

The open-circuit voltages on each short electrical 
dipole (or input voltages at each receiver channel) 
can be expressed as: 

Vi = hiE ( 1 ) 

where hi is an effective antenna length vector, E 
is the electric wave vector of the incoming wave and 
i denotes the antenna, i.e. i = X,Y, Z for DF2 mode 
or i = AY, Z and i = A, ZY for DF1 mode. Each h 
vector is represented by the following spherical coor¬ 
dinates: the antenna length hi and the colatitude Q 
and the azimuth £, angles. We chose the same coor¬ 
dinate system, as was used by Oswald et al. [2009] for 
rheometry and wire-grid calibrations of the SWAVES 
antenna. In particular, the colatitudes Cx,y,z of each 
SWAVES effective antenna vectors as well as the co¬ 
latitude of the AKR source are defined as an an¬ 
gle from the spacecraft axis A towards Y SC Z SC plane. 
The azimuth angles fx,Y,z and are counted from 
— Z sc towards Y sc . 

The sample covariance of the Vi in equation (1) 
which are actually measured by the receiver is lin¬ 
early related to the covariance of the incoming wave 
electric field (which basically forms the source polar¬ 
ization vector). Then, in the spacecraft coordinates 
the model-predicted analytic signals of the output 
voltages at each receiver can be written as [Ladreiter 
et al., 1995; Vogl et al., 2004; Cecconi and Zarka, 
2005]: 

5R«ViV/}) = 0.5Shihj[{l + Q^iClj + (2) 

H+ <>/!',) + (1 - 
Z{(ViV*)) = 0.5 ShihjVlQiVj -QjVi] 

where f= cos Q t j sin£ s — sin£j j:) - cos (j s cos(</> s — 
f'i.j ), 'll j_'j sin (jij sin(£ s ) ? C s and colat¬ 

itude and azimuth of the radio source in the space¬ 
craft coordinates; V):,j - analytical signals of the in¬ 
put voltages at each receiver; , Qj , fij - effective 
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length, colatitude and azimuth of electric antennas; 
S, Q,U,V - Stokes parameters of the incoming wave, 
asterisk denote complex conjugate and (...) is a time¬ 
averaging operation. 

During the spacecraft roll maneuvers SWAVES/ 
HFR1 was sequentially switched between DF1 T3’, 
DF1 ’31’ and DF2 modes for each measurement 
sweep. Therefore, the model-predicted signals are 
expressed the system of 4 equations in case of the 2 
antenna modes, DF1 T3’ and DF1 ’31’ or 9 equations 
for DF2 mode. The coefficients of these equations 
contain non-linear expressions of parameters describ¬ 
ing the the source direction (in case of point source 
approximation), Stokes parameters of the wave and 
the components of the antenna effective length vec¬ 
tors (colatitude, azimuth and length). 

With a distant observation at w 110 Re away 
from Earth the AKR radio sources can be assumed 
as a point source. During the roll maneuvers on 31 
Jan., 2007 STEREO-B was at high latitudes over 
the Northern hemisphere and therefore was able to 
observe only the Northern AKR sources which are 
mainly located in the evening sector of the magne¬ 
tosphere i.e. between ~ 18 : 00 — 22 : 00 hours of 
the magnetic local time (MLT) and ~ 65 — 75° of 
the invariant magnetic latitude(77anas^ et al. [2003]; 
Mutel et al. [2004]). In our calculations we assumed 
that the AKR source is located at 21:00 of MLT, 70° 
of invariant magnetic latitude and at the altitude of 
0.8 Re above the Earth surface (corresponds to the 
sources at ~ 300 kHz.) Therefore, the radio source 
coordinates in equation (2) are assumed as known 
parameters. Moreover, using the well known AKR 
polarization characteristic, i.e. that AKR is fully cir¬ 
cularly polarized (e.g. Panchenko et al. [2008]) with 
domination of the right-handed R-X mode (Hanasz 
et al. [2003]) as well as assuming that the polariza¬ 
tion state of the AKR does not change during the 
observation, we can suggest that the Stokes param¬ 
eters are Q = U = 0 and V = — 1. In addition, 
by assuming a perfect subtraction of the background 
contribution, one can see that the linear system of 
equations (2) is homogeneous and therefore it can 
be normalized to the unknown source intensity pa¬ 
rameter, and only the antenna effective length ratio 
retained. 

Due to the errors in the measurements we can¬ 
not expect an exact solution of the equations (2). 
Therefore the best solution in the least square sense 
is sought, i.e. the sum of the squares of the dif¬ 
ferences between the model predicted values (P mod ) 


and the observations (P ofes ) is minimized: 


M N / pobs (j. 

2 _ __2 I r m,n\ L n 

^ — Z_^ Z^ a m,n I Jobs 
m =1 n= 1 \ n 


Tmod 

± n 


2 


mi 43) 


where X = {( s , £s,hi,j,Q,j>€i,j>S, Q , U,V} are the 
input parameters of the model described by equa¬ 
tions (2)), index m counts the output receiver chan¬ 
nels (e.g. for DF2 mode: m = 1 denotes Ax, 
m = 2 corresponds to Ay, etc.), M is the to¬ 
tal number of the receiver output channels, index 
n denotes each data point, N is the total num¬ 
ber of data points in given time series and cr ra ,m 
represents uncertainties of each observation P° b ^. 
The normalization coefficients in equation(3) are : 

jobs, mod _ pobs, mod pobs, mod pobs, mod ^ pjp2 
mode, Iobs,mod = pobs^mU + pobs,mod for Dpl > 13 > 
mode and I° bs ’ mod = pobs, mod + p obp mod for Dpl 

’31’. 

Equation (3) can be used only in assumption that 
the errors of observations are uncorrelated. This 
implies that the variance-covariance matrix of ob¬ 
servations whose elements tr mj n are weights for the 
squared residuals is diagonal. In fact, when the an¬ 
tenna system receives the strongly polarized radio 
signal (e.g. AKR) one can expect cross correlation 
between uncertainties of the measurements and as 
a result the non-zero off-diagonal elements in the 
variance-covariance matrix. Therefore, in more gen¬ 
eral case methods of the variance-covariance matrix 
estimation described in Lecacheux [2011] should be 
used. In our calculations we used the observations 
with large signal/noise ratio and the expected devi¬ 
ations of the electric antenna vectors from orthogo¬ 
nality are not large. Therefore we assumed that the 
off-diagonal elements in variance-covariance matrix 
are very small and we can use the equation (3). 

The non-linear problem (3) can be solved by means 
of non-linear optimization methods. Most of these 
methods are based on iterative procedures: from a 
starting point xq the method produces a series of vec¬ 
tors Xi,X 2 , ...x n which converge to a '.true - the local 
minimizer for the function ~)c(x). The main prob¬ 
lem of such an algorithm is that the system of equa¬ 
tions can be ill-conditioned (numerically very close 
to singularity) in certain regions of the parameters 
space, with no unique solution. This is a well-known 
problem when the inversion is applied. In general, 
the ill-conditioned system of equations describes the 
situation when small fluctuations of the input data 
may result in a drastic increase of errors in the solu- 
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tion. The other problem of the nonlinear optimiza¬ 
tion method is the definition of the initial guess of the 
model parameters. Unfortunate choosing of the ini¬ 
tial guess can result in finding local minima of y 2 (%) 
which are much larger than its global minimum. 

3.3. Genetic algorithm (GA) 

Recently the stochastic global search techniques, 
such as genetic algorithms, became a powerful tool to 
solve inverse problems. The genetic algorithm (GA), 
the concept of which has been introduced by Holland 
[1975], is a stochastic technique based on the biolog¬ 
ical principle of survival of the fittest. The method 
represents each potential solution of the optimization 
problem as a set of chromosomes (genome). The al¬ 
gorithm begins from initialization of the initial pop¬ 
ulation of randomly selected ” individuals” (potential 
solutions). Each ’’individual” contains a set of chro¬ 
mosomes which represent the free parameters of the 
optimization problem. Applying the genetic opera¬ 
tors such as mutation and crossover the fittest ” indi¬ 
viduals” (potential solutions with the best approxi¬ 
mation of the problem) are determined on each iter¬ 
ation step. These fittest ’’individuals” have a greater 
chance to leave most offspring forming a next popu¬ 
lation to which the mutation and crossover operators 
are applied again. In this way, a successive approxi¬ 
mation of the solution of the optimization problem is 
achieved improving gradually. The application of the 
GA in astronomy and astrophysics is well reviewed 
by Charbonneau [1995]. 

The main advantage of the GA is that this algo¬ 
rithm doesn’t require a good initial guess of the solu¬ 
tion and therefore can operate in the whole space of 
the free parameters of the model. Therefore, unlike 
’’classical” nonlinear optimization methods, the GA 
is significantly less sensitive to misleading local min¬ 
ima and in most cases the global extremum of the 
function can be found. Genetic algorithms can also 
be adopted to solve ill-conditioned inverse problems 
[e.g. Mera et al., 2004], These advantages are coun¬ 
tered by fact that the GA requires higher computing 
power. Moreover, since GA is a stochastic method 
which uses random search, the resulting extremum of 
the optimization task is an extremum only in prob¬ 
abilistic sense. In other words, the obtained results 
with some probability are only approximations of the 
’’true” solution. The other source of difficulties in 
GA is the estimation of the solution errors. Gener¬ 


ally, in order to obtain the resulting uncertainties the 
stochastic techniques, e.g. ’’quick and dirty” Monte 
Carlo, are used [Press et al., 1992, pp. 689-699]. 

Charbonneau [1995] presents the examples of 
solving optimization problems by means of the 
GA as well as provides an open source opti¬ 
mization code PIKAI which maximizes a function 
f(xi,X 2 ,X 3 , ..., x n ) in the n-dimensional parameter 
space. PIKAI has the following major controlling 
parameters which determine the performance of the 
algorithm: number of generations N g (or number of 
iterations), size of the populations N p in each gen¬ 
eration, mutation rate p m and crossover rate p c . N g 
and N p are most critical parameters which determine 
the performance of the algorithm. The iterative GA 
can be terminated either after a given number of it¬ 
erations or when some convergence criteria are sat¬ 
isfied, e. g. improvement of the best solutions (the 
fitness level of the best individual) falls below a given 
threshold. In our case the size of the populations was 
N p = 100 and GA was terminated after N g = 300 
iterations. These numbers have been chosen as a rea¬ 
sonable compromise between computation time and 
accuracy of y 2 minimization. 

4. Results and error analysis 

As was mentioned in the previous section the GA 
can operate in the whole space of the free parameters 
of the model and does not require the initial guess 
of the parameters. Nevertheless, in order to reduce 
the computing time, we limited the space of possi¬ 
ble parameters in the following way: 1) tilt angles of 
each effective antenna vector h relative to physical 
antenna roods (see tables 1, 2, 3 columns ’’Mechani¬ 
cal” ) are limited to be 30°, and 2) the effective length 
of each antenna \h\ does not exceed ±50% offset from 
the length of the short dipole antenna (L/ 2). 

The red lines in Figure 2 show the model predicted 
values for DF2 mode fitted to the AKR observations 
using GA. As seen, the model fit the observations 
well. The resulting effective antenna lengths for dif¬ 
ferent antenna configuration (optimized model pa¬ 
rameters) are listed in tables 1, 2 and 3. 

One important point is that since measurements 
are contaminated by instrumental systematic errors 
as well as a random noise the best fitted set of pa¬ 
rameters is not a unique realization of the ” true” pa¬ 
rameters. Therefore, the errors of the fitted param¬ 
eters have to be carefully estimated. In a nonlinear 
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STEREO-B/WAVES, DF2 mode, 01/31/2007 04:30 - 15:00 UT 
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Figure 2. Modeled analytic signals (red lines) fitted to the measured output voltages at each 
channel of the STEREO-B HFR1 receiver working in Direction Finding (DF2) mode (blue lines). 
The modeled and measured time profiles are normalized by I = Ax + Ay + Az- In DF2 mode 
all three antennas (X, Y, Z) are operated. The best fitted parameters which define the effective 
antenna vectors are tabulated. 


least-squares methods (e.g. conjugate gradient-type 
techniques) the errors can be obtained from inverse 
Hessian matrix of second derivatives evaluated at the 
best solution. Since GA does not use the derivatives 
the best way is to use a Monte Carlo simulation in 
which a large number of synthetic data sets are gener¬ 
ated from the best fitted model, adding random noise 
in accordance with the measurement errors. The GA 
algorithm is applied for each of these sets and finally, 
the standard errors of the obtained estimators are de¬ 
termined. 

In order to estimate the fit errors, we used a so 
called ’’Quick and dirty” Monte Carlo or Bootstrap 
method (e.g. Efron and Tibshirani [1993]) - a power¬ 
ful technique which can be applied even without in¬ 
formation about a true underlying error distribution 
of the observations. The bootstrap method generates 
the synthetic data sets by randomly resampling with 
replacing of approx 37% of the original data (see for 
more details [Press et al., 1992, pp. 689-699]). We 


synthesized M=500 bootstrapped data sets using the 
measured temporal variations of the AKR intensity 
and then applied the fitting procedure based on GA 
to each of these synthetic data sets (each synthetic 
set has the same number of points as the measured 
AKR intensity profile). It provides 500 sets of pa¬ 
rameters Xi, ...,X m each of which define the effec¬ 
tive antenna length vectors. Figure 3 shows the his¬ 
togram of the bootstrapped parameters. Then, since 
the distribution of the bootstrapped parameters Xj 
is approximately normal (as seen in Figure 3), the 
widths of the distributions give us the standard de¬ 
viations a of the estimated parameters. 

The same calculations have been performed also 
for the SWAVES antennas working in DF1 mode and 
the final results with one sigma errors are listed in 
tables 1, 2 and 3. 

5. Discussion and summary 
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Figure 3. Bootstrap error analysis of best fitted parameters X = 

{hx,Cx,^x',hY/hx,CY,^Y',hz/hx,Cz,^z} which define the effective vectors of the SWAVES 
antennas for DF2 mode. The histograms show the distributions of the fitted parameters obtained 
by application of the genetic optimization algorithm to M=500 synthetic data sets generated by 
means of the bootstrap method. The standard deviations of the X elements are the same as the 
standard deviations a of the bootstrap histograms. More details are given in the penultimate 
paragraph of the section 4. 


Table 1. Results of in-flight calibration of the SWAVES antennas in direction finding mode DF2 
(monopoles X , Y, and Z). The column headed ’’Mechanical” contains relative lengths and directions of 
the antenna rods. The column headed ’’Effective” contains the components of the effective antenna lengths 
and directions with one sigma errors. £ and £ are defined in the STEREO/WAVES spherical coordinates, 
described in section 3.2, and 9 and p are colatitude and azimuths of the antenna vectors in the spherical 
coordinates related to the spacecraft main axis, i.e. 9 is measured from Z spacecraft axis, and ip is azimuths 
angle from X spacecraft axis. 




Mechanical 


Effective 


Ant. 

h/hx 

c 

C 

h/h x 

c 

€ 

X 

1.0 

125.26° 

-120.0° 

1.00 

120.8° ±0.5° 

-135.6° ±0.8° 

Y 

1.0 

125.26° 

120.0° 

1.27 ±0.01 

116.7° ±0.3° 

124.7° ± 0.7° 

Z 

1.0 

125.26° 

O 

b 

o 

0.78 ± 0.01 

124.9° ±0.5° 

13.3° ±0.8° 

Ant. 

h/h x 

9 


h/h x 

9 


X 

1.0 

65.9° 

230.8° 

1.00 

52.1° ±0.7° 

229.6° ± 0.6° 

Y 

1.0 

65.9° 

129.2° 

1.27 ±0.01 

59.4° ±0.6° 

121.5° ± 0.4° 

Z 

1.0 

144.7° 

180.0° 

0.78 ±0.01 

143.0° ± 0.5° 

161.8° ± 1.0° 


The final set of antenna parameters of STEREO/ 
WAVES obtained by the in-flight calibration (for di¬ 
rection finding mode DF2) are compared with rheom- 
etry and computer simulations previously performed 


in Graz and reported in Macher et al. [2007] and 
Oswald et al. [2009]. Table 4 summarizes the ele¬ 
ments of the antenna effective length vectors using 
the different methods. For each method we have 
also checked the goodness of the fit by calculating 
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Table 2. Results of in-flight calibration of the SWAVES antennas in mode DF1 ”13” (monopole X - dipole ZY). 




Mechanical (’13’ mode) 

Effective (’13’ mode) 

Ant. 

h/h x 

c 

1 

h/hx 

C 

€ 

A 

1.0 

125.26° 

-120.0° 

1.00 

120.4° ±0.8° 

-136.3° ± 1.1° 

ZY 

1.414 

90.0° 

150.0° 

1.32 ± 0.01 

92.4° ± 1.0° 

146.7° ± 1.0° 

Ant. 

h/hx 

e 

f 

h/h x 

e 


X 

1.0 

65.9° 

230.8° 

1.00 

51.4° ±0.9° 

229.7° ± 0.9° 

ZY 

1.414 

30.0° 

90.0° 

1.32 ± 0.01 

33.4° ± 1.0° 

94.4° ± 1.8° 


Table 3. Results of in-flight calibration of the SWAVES antennas in dipole mode DF1 ”31” ( monopole Z - dipole XY ). 


Mechanical (’31’ mode) Effective (’31’ mode) 


Ant. 

h/hz 



Whz 



XY 

1.414 

90.0° 

90.0° 

1.76 ±0.01 

90.6° ±0.7° 

88.7° ± 0.8° 

Z 

1.0 

125.26° 

O 

o 

o 

1.00 

124.3° ± 0.7° 

13.8° ± 1.0° 

Ant. 

h/hx 

e 

v 

h/h x 

e 

<P 

XY 

1.0 

90.0° 

90.0° 

1.00 

91.3° ±0.8° 

90.6° ± 0.7° 

ZY 

1.414 

144.7° 

180.0° 

1.33 ±0.01 

143.3° ± 0.7° 

160.7° ± 1.3° 


the reduced chi-squared: £ = % 2 /{N * M — u — 1), 
where y 2 is defined by equation (3), N is the total 
number of observation sets, M is the number of out¬ 
put channels, and v the number of fitted parameters 
of the model. £ describes the discrepancy between 
the observed signals of the AKR and the model pre¬ 
dicted values, as calculated by equation (3) with the 
effective antenna length vectors obtained by the re¬ 
spective method. As can be seen from the table, 
the smallest £ corresponds to the in-flight results, al¬ 
though the values of the goodness of the fit (£) are 
very similar for all methods. 

Also, the results of the experimental rheometry 
and numerical computer simulations well agree with 
the in-flight calibrations analysis within the inherent 
inaccuracy of the methods. In particular, the angu¬ 
lar components of the antenna effective length vec¬ 
tors obtained by the different methods are within 3cr 
errors of each other. This confirms the high reliabil¬ 
ity of the rheometry and computer modeling for the 
calibration of the antenna system. For the effective 
lengths we find very small but significant differences 
between the present results and those obtained by 
former authors ( Macher et al. [2007]; Oswald et al. 
[2009]; Zaslavsky et al. [2011]). The most plausible 
explanation is that the base capacitances associated 
with the antennas X , Y and Z are not exactly the 
same. Since base capacitances shorten the effective 
length of an antenna, the results given in Table 4 
indicate Cy/Cx < 1 and Cz/Cx > 1 (the latter 


being confirmed by Zaslavsky et al. [2011], who have 
determined hx and hz, but not hy). 

As was mentioned in section 3.1 only STEREO-B 
observed the AKR during the roll maneuvers and, 
therefore, the in-flight calibration was applied only 
for the SWAVES antenna system onboard STEREO- 
B. Nevertheless, since the two STEREO spacecraft 
are almost identical we expect that the antenna ef¬ 
fective length vectors of the WAVES experiment on 
STEREO-A are very close to STEREO-B, deter¬ 
mined by the in-flight calibration. This assumption 
was confirmed by rheometry and computer simula¬ 
tions. However, Zaslavsky et al. [2011] found slightly 
different results for STEREO-A and STEREO-B by 
analysis of the effective antenna lengths using the 
galactic background. 

The proper calibration of the reception properties 
of the spaceborne antenna system is one of the major 
data processing tasks yielding an accurate evaluation 
of the radio data, which is very important especially 
for the direction finding of electromagnetic waves. 
The in-flight calibration approach, presented in this 
paper, gives us reliable results because this method 
uses the data acquired by the ” real” antenna-receiver 
configuration which are difficult to estimate or mea¬ 
sure during the spacecraft ground-based calibration 
procedures. For example, the accurate determina¬ 
tion of the stray and base capacitances of the an¬ 
tennas plays an important role in achieving realistic 
results by rheometry and computer simulation meth¬ 
ods [ Macher et al., 2007]. 
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Table 4. Results of the STEREO-B/WAVES antenna calibration by in-flight calibration, computer simu¬ 
lation, and rheometry (computer simulation with CONCEPT and for ” loaded” feeds and rheometry results 
from Oswald et al. [2009]). Column ’’Mechanical” contains the values which describe the relative length and 
direction of the antenna rods. The values in the parentheses denote the standard deviation of the obtained 
values. £ is the reduced chi-squared statistic which is to check the goodness of the fit. 


Ant. 


In-flight 

Computer simulation 

Rheometry 

Mechanical 


hx/h x 

1.00 

1.00 

1.00 

1.00 

X 

4, deg 

120.8° (0.5) 

119.9° 

121.3° 

125.3° 


4, deg 

-135.6° (0.8) 

-135.3° 

-135.4° 

-120.0° 


hy/hx 

1.27(0.01) 

1.21 

1.26 

1.00 

Y 

4,, deg 

116.7°(0.3) 

114.4° 

115.1° 

125.3° 


4, deg 

124.7°(0.7) 

127.3° 

126.8° 

120.0° 


h-z/h x 

0.78(0.01) 

0.81 

0.84 

1.00 

Z 

4, deg 

124.9° (0.5) 

124.7° 

125.3° 

125.3° 


4, deg 

13.3°(0.8) 

15.5° 

16.4° 

0.0° 

4 


1.08 

1.26 

1.25 

11.3 


Since the AKR is observed in the frequency 
range below 1 MHz, we can only study the recep¬ 
tion properties of the antenna in electrically short 
dipole approximation when the radio wavelength is 
much greater than the size of the antenna system 
- i.e. Lantenna < >^wave- In this frequency range 
(also called quasi-static frequency range) the effective 
length vectors are real and constant (independent of 
direction). The effective antenna length vectors on 
higher frequencies will be significantly different, espe¬ 
cially for frequencies close to the antenna resonances 
[Oswald et al., 2009]. 
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